Quantum wire networks with local Z2 symmetry 
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For a large class of networks made of connected loops, in the presence of an external magnetic 
field of half flux quantum per loop, we show the existence of a large local symmetry group, generated 
by simultaneous flips of the electronic current in all the loops adjacent to a given node. Using an 
ultra-localized single particle basis adapted to this local Z2 symmetry, we show that it is preserved 
by a large class of interaction potentials. As a main physical consequence, the only allowed tunneling 
processes in such networks are induced by electron-electron interactions and involve a simultaneous 
hop of two electrons. Using a mean-field picture and then a more systematic renormalization-group 
treatment, we show that these pair hopping processes do not generate a superconducting instability, 
but they destroy the Luttinger liquid behavior in the links, giving rise at low energy to a strongly 
correlated spin-density-wave state. 

PACS numbers: 71.10.Li,71.10.Pm,73.21.Hb 



I. INTRODUCTION 

For the past two decades, transport properties of quan- 
tum wires have received a lot of attentionii^. Besides 
metallic systems, two dimensional electron gases induced 
in GaAs/AlGaAs heterostructures have displayed a rich 
variety of quantum effects such as Aharonov-Bohm re- 
sistance oscillations in a ring geomctrji^, and persistent 
current9^. More recently, coherent Aharonov-Bohm os- 
cillations have been measured in ballistic arrays with the 
dice lattice geometrjii, in agreement with the predictions 
of simple models for non-interacting electrons^. The ex- 
perimental visibility of conductance oscillations with the 
full flux quantum $0 = /i/e periodicity in large struc- 
tures with a few thousand loops has been attributed to 
the presence of Aharonov-Bohm cages which localize the 
motion of electrons when the external flux per elementary 
loop is half-integer in units of <i>o. 

After these first investigations, a natural issue is to un- 
derstand interaction effects on this geometrical localiza- 
tion phenomenon. For two electrons, it was shown that 
interaction induces some correlated hopping processes, 
at the origin of some delocalized states in which the two 
electrons always remain close from each otheri2ii£. Re- 
markably, this trend appears for any sign of interaction. 
Of course, for a repulsive interaction, these extended 
states lie above the ground-state sector, obtained when 
the two electrons are individually localized in distinct re- 
mote cages. But this raises the question of what happens 
at finite density, where it is no longer possible to put all 
conduction electrons in separate cages. Two electrons 
initially in the same cage will be able to move together, 
so this may generate some kind of metallic state. 

An interacting system with a finite particle density was 
first studied in its bosonic version in a one-dimensional 
geometry showing Aharonov-Bohm cages, namely a chain 
of rhombiii. This model has been adapted into a real- 
istic description of a Josephson junction array with the 



same geometry recentljii^. There, it was found that no 
single boson condensate can exist, but instead a 2-boson 
condensate is stabilized at low energy. This peculiar phe- 
nomenon is connected to the single particle localization 
in Aharonov-Bohm cages, but in the presence of interac- 
tion, it requires a stronger property, identified as a local 
symmetry based on the Z2 group. This symmetry im- 
plies that the parity of the total particle number within 
any Aharonov-Bohm cage is conserved. If this symmetry 
remains unbroken, it prevents the formation of any single 
boson condensate. A two-dimensional version exhibiting 
this local Z2 symmetry was also analyzedi^, in the con- 
text of topologically protected quantum memories. The 
crucial role of a local symmetry is well illustrated by the 
example of the dice lattice, which sustains Aharonov- 
Bohm cages, but no local symmetry in the presence of 
interactions. There, it was found that the external mag- 
netic field corresponding to <I>o/2 per loop imposes a vor- 
tex lattice which allows for some zero-energy line defects, 
but the single particle condensate exhibits a finite phase 
stiffnesai^ii^. 

The fermionic case at finite density is the subject of 
the present article. Motivated by ongoing experimentsii 
on networks designed on a 2D electron system at a 
GaAs/GaAlAs interface, we shall study quantum wire 
networks, which geometry is compatible with a local Z2 
symmetry in the presence of a magnetic field of half flux 
quantum per elementary loop. Such networks will be 
called Z2 networks, and two examples are shown on Fig.^ 
below. In section |nl we define precisely this class of net- 
works, construct the unitary operators associated to the 
local Z2 symmetry group, and consider their irreducible 
representations in the Hilbert space of a single particle. 
We show that all the common eigenstates of all the lo- 
cal generators of this large Abelian group are strongly 
localized within Aharonov-Bohm cages involving only a 
finite number of elementary loops. Section ITTll first dis- 
cusses the matrix elements of a two-body interaction in 
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FIG. 1: This picture presents two examples of Z2 networks: 
a chain of loops (upper part), and a square lattice of loops 
(lower part). All these networks are obtained from any planar 
lattice after replacing each link by a loop. To get a local 
Z2 symmetry, each loop is required to be perfectly invariant 
under the reflection through the long axis defined by its two 
adjacent nodes. 

this localized basis adapted to the Z2 symmetry. This 
shows that a large class of pair interaction potentials (in- 
cluding the point-like Hubbard interaction) is compat- 
ible with this symmetry. We then show that, at least 
in a mean-field approach, the particle pair hopping pro- 
cesses generated by a repulsive interaction do not lead to 
superconductivity. Starting from simple order of magni- 
tude estimates adapted to semi-conductor quantum wire 
networks, we argue that a more refined renormalization 
group approach (RG) is required, given the strength of 
the Coulomb interaction for typical electronic densities. 
Finally, section IfVI proposes a simplified Z2 symmetric 
model of coupled Aharonov-Bohm cages which is ana- 
lyzed with the RG method using the one-loop approxi- 
mation. This confirms the irrelevance of the electron pair 
hopping processes in the low-energy limit, and establishes 
a tendency for the system to reach a strong coupling fixed 
point of a spin density wave character. So the combina- 
tion of local Z2 symmetry and electron-electron repulsion 
destroys the Luttinger liquid in the quantum wires which 
constitute the network, leaving room only for a strongly 
correlated liquid with mostly gapped excitations. Finally, 
a brief appendix shows how this picture fits with a simple 
mean-field description. 



II. LOCAL Z2 SYMMETRY FOR A CLASS OF 
QUANTUM NETWORKS 

A. Basic construction 

In this paper, we shall consider networks of quantum 
wires (called Z2 networks), for which the elementary links 



connecting two nearest-neighbor nodes are doubled. In 
other words, such networks may be obtained from any 
planar graph, after replacing links by loops. Examples 
of this construction are shown on Fig. ^ The key as- 
sumption we shall make is that all these loops exhibit a 
perfect refiection symmetry around the axis joining the 
two nodes connected by the loop. This requirement im- 
poses some rather severe constraints on the fabrication 
of such structures, in terms of the network geometry, 
and of the control of various types of disorder. However, 
as discussed in the Introduction, experimental observa- 
tion of the Aharonov-Bohm cage effect in GaAs quantum 
wire networks^ suggests that present technology is able 
to meet these constraints. On a more theoretical side, 
numerical investigations on disordered tight-binding ver- 
sions of such networks have shown that the effects we 
shall discuss survive after a small but sizeable disorder 
is introducediS,. At a purely geometric level, these struc- 
tures exhibit, besides their usual crystalline symmetry 
groups, a large symmetry generated by local flips of ele- 
mentary loops around their central axis. We are mostly 
addressing here the quantum-mechanical effects of such a 
symmetry in the presence of a constant and uniform ex- 
ternal magnetic field, corresponding to half-fiux quantum 
per loop. 

As usual, quantum mechanics is affected by the fact 
that the vector potential used to describe the magnetic 
field is not invariant under these local operations. In 
the case of the translational symmetry, this leads to the 
non-commutation of magnetic translations and to Lan- 
dau level quantization. In the present case, the striking 
consequence is the appearance of a dramatic localization 
phenomenon, where all single particle energy eigenstates 
are confined inside a finite set of loops^. 

To illustrate the idea, it is useful to consider first a 
single loop, with strictly one-dimensional motion along a 
coordinate x. If / denotes the perimeter of this loop, we 
have to impose the periodicity condition: ipil) = '0(0). 
The stationary Schrodinger equation reads: 

Eip{x) = -^(p + eAfi}j{x) = ^{-ifi-^ + eA)^V(a;) 
2m 2m ax 

where A(x) is the projection of the vector potential 
along the wire at point x. Using the gauge transforma- 
tion: -0(a;) — cxp{—i^ A{x')dx')xix), xi^) obeys the 
Schrodinger equation for a free particle, subjected to the 
boundary condition: %(/) = exp(i27r$/<I>o)x(0)! where 
$ = I A{x')dx' is the magnetic fiux through the loop 
and <I>o = h/e is the fiux quantum. The single particle 
energy spectrum is then: En = ^^^^-(ri 4- $/$o)^7 where 
n is any integer. For generic values of ^/^o- Sill energy 
eigenvalues are non-degenerate, but all levels are doubly 
degenerate when 2$/(f>o is an integer (with the exception 
of the ground-state for integer 4>/<I>o). For integer <i>/<I>o, 
the flux through the loop may be eliminated by a gauge 
transformation, and the degeneracy is clearly attributed 
to the reflection symmetry, which changes the wave- 
number of an eigenstate into its opposite. In the case 
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of half- integer $/$0i this geometrical symmetry has to 
be combined with a gauge transformation so that ip^x) is 
changed into {Tij}){x) = tl''{x) ~ cxp{if {x))ip{~x) , with 

f{x) = f (/q"'' A{x')dx' - A{x')dx'^ . This transfor- 
mation sends an eigenstate with a persistent current into 
an eigenstate with the reversed current. We also notice 
that ^/y{0) = -0(0) whereas ^p'{l/2) = -iil/2). This has 
some important consequences if we wish to generalize 
these reflection operations to the elementary loops con- 
necting nodes in Z2 networks. 

Let us first consider a one dimensional chain of loops, 
as shown on the upper part of Fig. ^ Any link of this 
chain can be represented by a simple line with plus and 
minus signs on its extremities, see Fig. |21 reminding us 
the sign of the transformed wavefunction. Suppose now 
we try to construct a symmetry associated to the reflec- 
tion for just one loop. The simplest way is to combine 
as before reflection and gauge transformation, but we en- 
counter a major problem: this transformation will be no 
longer local, because we are obliged to change globally 
the sign of the part of the wave-function lying on one 
side of the transformed loop. For an open chain, this 
non local aspect of the transformation attached to a sin- 
gle loop is not too disturbing, since it modifies physical 
observables such as electronic current only locally. But 
in the case of periodic boundary conditions, it has the 
consequence that such a transformation simply becomes 
impossible. In fact, one can easily avoid this difficulty by 
considering two transformations on adjacent loops, see 
Fig. 13 The generalization for any Z2 lattice is clear, we 
need to join together negative extremities of all the loops 
arriving at a given node. Because of these phase-factors 
arising from gauge transformations, the local symmetry 
generators are therefore associated to lattice sites and not 
to the links as one would infer from purely geometrical 
considerations. Note that in lattice gauge theories, local 
generators of gauge transformations are also attached to 
lattice sites, which establishes a close connection between 
Z2 networks and lattice gauge theories based on the Z2 
group. This viewpoint has been emphasized in the con- 
text of Josephson junction arraysi^ii^. 

In case of an integer number of flux quanta per loop 
($ = n$o) , we have now ^'(O) = -0(0) and 0'(V2) = 
ip{l/2), so that each link carries plus signs on both ends. 
In this case, it is possible to construct local Z2 symmetry 
generators associated to lattice links instead of nodes. 
But then, because the single particle Hamiltonian has 
mostly a dispersive spectrum, simultaneous eigcnstates of 
the Hamiltonian and these local generators are extended, 
instead of being localized as in the half-flux case. 



B. Constraints on scattering matrices at the nodes 

Let us now describe more precisely the local transfor- 
mation attached to a given site i on a Z2 network. This 
site belongs to n elementary loops, so it is the intersec- 




FIG. 2: This picture shows different realizations of a local 
Z2 symmetry, successively for a single loop (upper part), and 
then around sites present in a ID geometry (chain of loops, 
lower left) or a 2D geometry (square lattice of nodes con- 
nected by loops, lower right). Signs plus or minus keep track 
of the relative sign between initial and transformed wavefunc- 
tions ■(/) and T-i/) = -(/''• 

tion of 2n links. Since these links are naturally paired, 
we shall label them by a compound index (ija), where 
J refers to any of the n nearest neighbors of site i, and 
a £ {1,2}. On each of these 2n links, we choose the 
origin of an x coordinate at site i, and the projection of 
the local vector potential along this link is denoted by 
Aija{x). As already discussed, it is convenient to write 
the single-particle wave-function as: 

■ipijaix) = exp ^-i^ J Aija{x')dx'^ Xija{x) 

along the link (ija), where an eigenstate with energy 

1-2 I 2 

E = has the form: 

\ija{x^ — Cijae ~\- Dijae 

Here, Cija (resp. Dija) is the amplitude of an incom- 
ing (resp. outgoing) wave with respect to site i. Since 
Xijai^) is locally the wave-function of a free particle with 
no magnetic field, we may define in a gauge-invariant way 
a scattering matrix at site i by: 

-^y" ^ X! ^jl,j'a''^i3'a' (1) 
j'a' 

which expresses outgoing amplitudes as a function of in- 
coming ones. Using the notation 1 = 2 and 2 = 1, the 
local Z2 generators commutes at any site with the single 
particle Hamiltonian if and only if: 

c(i) ^ c(i) ^ oii) / ■ / (0) 

jaj'a' ja.j'a' ja,j'a' \J / J / V / 

ja,ja' ja,ja' \ ) 

where i denotes any site, and j, j' are arbitrary nearest 
neighbors of i. In other words, there should be a perfect 
symmetry between paired links ija and ija. The proof 
of this rather intuitive statement will be given in the 
next paragraph. Note that in real systems, nodes always 
have a finite area, and this constraint may be difficult to 
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implement. For instance, as shown on Fig. in a chain 
of loops, it seems very likely that the simplest design for 
a node docs not satisfy this symmetry. A suggestion for 
improving this local geometry is sketched on the bottom 
of Fig. 121 Generalizations to multichannel quantum wires 
are possible, if the underlying inversion symmetry is also 
performed on the transverse channels. But we have to 
restrict ourselves to a regime of weak magnetic fields so 
that we may neglect the field effects on orbital motion 
within a wire. 






FIG. 3: This picture illustrates the requirements on the scat- 
tering matrix at a node present in a chain of loops (upper 
part), in order to impose a local Z2 symmetry. The most 
direct design (a), is likely to be incompatible with these con- 
straints, since intuitively, the amplitude for an electron to go 
from A to B is different from the amplitude to go from A 
to C. In a slightly modified design (b), these amplitudes are 
more likely to be equal, and this is more appropriate for an 
experimental implementation of the local Z2 symmetry. 

Let us now check that conditions ©j® are all we 
need. The transformation attached to site i now reads: 

i''^JC.{^) - (4) 

- exp (^<^a^ j (.A^ij2(.x') - Ayi(a;'))(ix'^ Ajs{x) 

where cti = 1 and <T2 = ^1- Using the locally gauge- 
invariant description in terms of Xijaix), this trans- 



lates into: x'ijai^) — ~Xijaix), which is equivalent to 
C'lj^ ~ —Cijsn and -D^jQ, = —Dija- This operation leaves 
the S matrix at site i invariant if and only if: 



S, 



jaj'a' 



(5) 



Note that this is valid for any j and j' which are nearest 
neighbors of i. When j = j' , we recover condition I®, 
and when j 7^ j', it is implied by condition But 
we should also check that the transformation attached to 
site i also preserves the S matrix at all nearest-neighbor 
sites J of i. Along the link ija, we have to perform a 
change of origin, permuting i and j, and replacing x by 
y = lij —X, where Uj is the common length of the two links 
{ijl) and {ij2). Using the fact that the flux through the 
loop enclosed by these two links is half a flux quantum 
modulo $0, Eq. Q becomes: 



(6) 



exp 



3i2{y') - Aja{y'))dy' Vjm(y) 



which implies x'jtai^) = Xjiaix) 



Djia ■ This yields S. 



or: q„ 



and 



S: 



that is Eq. (0) 



after a permutation of i and j. Now, for any site j', 
nearest neighbor of j and distinct from i, the trans- 
formation based at site i does not affect the pair of 
links (jj'a), so C^,^ 



implies 



Cjj'a. and D'..,^ 



Djjia, which 



and S^i,^^^ = S^i,^^^,thatis 



Eq. |(5J) after a permutation of i and j 



C. Irreducible representations of local Z2 symmetry 

Let us denote by Ui the unitary transformation at- 
tached at site i, defined by Eq. Q). Clearly, Uf = 1. so 
the eigenvalues of Ui are ±1. Furthermore, it is simple 
to check that Ui and Uj commute. The only non-trivial 
case arises when i and j are nearest neighbors. In fact, 
comparing Eqs. |0J and (j^J) we see that the actions of 
Ui and Uj on the loop between i and j differ only by 
a sign. The mutual commutation of these operators en- 
sures that they can be simultaneously diagonalized, so 
that Ui\'^) = fil^"), with = ±1. For single particle 
states, the possible collections of eigenvalues {e^} are very 
limited. To identify them, we need to introduce the no- 
tion of an Aharonov-Bohm cage. For any lattice site i, 
the cage attached to i is composed of all the elementary 
loops which contain i. We have now the following prop- 
erties: 

Property (a): If i and j are nearest-neighbor sites, and 
if Uil"^) ~ Uj\'i>), then the wave- function associated to 
I'?) vanishes on the loop joining i and j. 
This is an immediate consequence of the fact that along 
this loop, Ui\'^} and J7j|^) are described by opposite 
wave- functions, as already noted (compare Eqs. Q and 



®) 
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Property (b): If = — l^*), the wave-function ip at- 

tached to I ^P) vanishes everywhere outside the Aharonov- 
Bohm cage attached to i. 

This follows directly from the fact that the transforma- 
tion Ui does not modify the wave-function outside the 
Aharonov-Bohm cage attached to i. 
Main Property: The only possible collections of si- 
multaneous eigenvalues {ei} are such that exactly one 
of them is equal to —1, all the others being equal to 
1. If i is the only site with = —1, the corre- 
sponding wave-function vanishes everywhere outside the 
Aharonov-Bohm cage attached to i. 
Because of property (a), if all e^'s are equal to 1, |^) van- 
ishes. So at least one eigenvalue is negative. Suppose we 
have two of them: et = ej = —1 with i ^ j. Because of 
property (b). the wave-function ip vanishes everywhere 
outside the intersection of the cages centered at i and j. 
If this intersection is empty, then j^") = 0. If not, then 
i and j are nearest neighbors. But then, property (a) 
implies that ip vanishes also on the loop joining i and 
j, so 1^) ~ 0. The only possibility is therefore to have 
ei ~ —1 and ej = 1 for any j ^ i. From property (b), 
the wave-function is confined inside the Aharonov-Bohm 
cage centered around site i. 

When the S matrix of each node satisfies conditions Q 
and lj2J) for any wave- vector k, we may diagonalize 
the single-particle Hamiltonian in a localized basis of 
Aharonov-Bohm cage wave- functions. Since each of these 
cages has only a finite total wire length, the energy spec- 
trum is always discrete. These cages behave as artificial 
atoms, at least as long as electrons are not interacting. 
The interesting feature of such a system is that two cages 
may overlap, which allows for non-trivial interaction ef- 
fects. 

All these results (for $ = $o/2) have a simple physical 
interpretation in the light of the Aharonov-Bohm effect. 
Suppose we have an electron situated at a node i. If j is 
a nearest-neighbor of i, there are two possible equivalent 
ways for this electron to reach a nearby loop {j, k), using 
either side of loop But the associated probability 

amplitudes have opposite signs (as exp(i27r$/$o) = ^1) 
and will compensate each other. This is why we obtain 
localized solutions which form a complete basis of the 
Hilbert space for non-interacting electrons. 

In case of integer number of flux quanta per loop, the 
functions within a given irreducible representation are 
not necessarily local, so the previous conclusions do not 
apply. It happens that some localized solutions still exist, 
but they form only "half" of a complete basis of states. 
This analysis could be made using the Kottos and Smi- 
lansky formalismSS, but to save space and since we are 
mostly interested in half-flux quantum per loop, we will 
not give further detail on the integer flux case in this 
paper. 



III. INTERACTION EFFECTS ON 
AHARONOV-BOHM CAGES 

A. Hamiltonian for interacting electrons in an 
Aharonov-Bohm cage basis 

In the case of interacting electrons, this simple phys- 
ical picture has to be deeply modifled, because the two 
sides of a given loop are no longer equivalent if an elec- 
tron "sees" on one of them another electron and interacts 
with it. So the probability for one electron to hop from 
an Aharonov-Bohm cage to a nearby one no longer van- 
ishes. However, this is an interaction-induced effective 
hopping, and it is by no means obvious that it could 
lead to a coherent metallic system in the usual sense. By 
contrast to lattices such as the dice lattice which also 
exhibit the Aharonov-Bohm cage phenomenon, we will 
now show that on Z2 lattices, and for a large class of 
electron-electron interactions, the local Z2 symmetry is 
still present. This puts some rather severe constraints 
on the amount of delocalization that interactions may 
induce for a finite density of electrons. Because of the 
local Z2 symmetry, only pairs of electrons are able to 
hop from one cage to an adjacent one. So on intuitive 
grounds, a true conducting state is expected only for at- 
tractive interactions, where it is also a superconductor. 
For repulsive interactions, we shall see that the system 
remains an incoherent metal. 

For the sake of simplicity, we shall consider here a local 
Hubbard repulsive interaction between electrons. Gener- 
alization to a large class of interaction potentials is given 
in Section IIII Bl below. This interaction may be written 
as: 

i^int f dx (a;)V'ljcx (^)'P^jc.\ i^) 

Jo 

(7) 

Here ipija-^ i^) is the second quantization electron creation 
operator at point x on the wire a = {1,2} joining nodes 
i and j, with spin projection "f along the quantization 
axis. Let us now introduce operators of creation and 
annihilation of an electron in a given Aharonov-Bohm 
cage state: cj^(a;) and Cj^(a;) respectively, r referring to 
the spin projection. They satisfy the following relations: 

[/.4(.t)[/, ^ -clix) (8) 
UjcU^W, = clix) i^j (9) 

which show that these operators c\^{x) and c^^{x) carry 
a non-trivial local Z2 charge. The precise correspondence 
between the initial operators i^lj^-^ix) and the cage op- 
erators c\.j.{x) is as follows: 

4^^,{x) = ^ (e-^-°(-'4(x) + e-^^-(^)c],(a;)) 

(10) 

where 6ija{x) = f /j AyQ.dr is a line integral oriented 
from node i to node j. The local Hubbard interaction in 
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cage basis takes the following form: 

4t (2;)ca {x)Ci^ {x) + 4^ (x)c]^ {x)Cj^ {x)c^^ {x) + [a] 

4t (^)S-i (^)'^»T + 4t (^)'^a (-^^S'T + W 
4t (^)^a (^)St + 4t (^')"^a (^)Si ^-^^^^t + M 

4t (2;)cji (^)st + 4t (^)^u (^)^*T 

Lines on the right-hand side represent respectively: an 
intra-cagc Hubbard interaction [a], an inter-cage Hub- 
bard interaction [b], a neighbor cage exchange process 
[c] , a pair propagation process [d] . 

In the case of a one-dimensional chain of loops, these var- 
ious processes may be visualized directly on the graph 
shown on Fig. 01 

A 

m 




FIG. 4: In the case of Hubbard like interacting electrons, 
the only possible cage transitions preserving local Z2 symme- 
try can be plotted on the following graph (thin lines), where 
(n, m) and (n ,m') indicate initial and final cage positions, 
corresponding to a term c^^^,c'^,CmC„. Dashed lines refer to 
possible non zero amplitude transitions, compatible with the 
local nature of two-electron interaction, but which are forbid- 
den since they violate electron number parity conservation 
resulting from the local Z2 symmetry. 

The main feature appearing in this new basis is that 
the parity of the total electron number in each Aharonov- 
Bohm cage is separately conserved. Because an electron 
in state c\^{x)\Q) carries a non-trivial Z2 charge, but a 
pair c\^{x)c\^{x)\Q) belongs to the identity representa- 
tion of the local Z2 group, the previous statement simply 
shows that the local Hubbard interaction commutes with 
the local Z2 symmetry. In Scction rill Bl we show that any 
pair potential interaction of the form F(r, r') is compat- 
ible with the local Z2 symmetry provided ^(r, r') does 
not change when r or r' is replaced by its mirror image 
with respect to the axis joining its two nearest nodes. In 
a previous work, a similar statement was reached for a 
rather special case, namely a chain of rhombi in a tight- 



binding descriptionii. The present formulation is much 
more general. 

What arc the consequences of this local symmetry on 
an interacting system? As ir&iS, let us first consider 
a pair of electrons, with the local Hubbard interaction. 
We have three possible situations. First, the two elec- 
trons may be located in two cages i and j which do not 
share common links, that is i and j are not nearest neigh- 
bors. Then they do not "see" each other and remain 
localized in their respective cages. A second possibility 
occurs when i and j are distinct but nearest neighbors. 
As term [c] in Eq. l(TT)l above shows, the two electrons 
can exchange their cages, but such state remains local- 
ized under time evolution. Finally, when i = 7, term [d] 
in Eq. indicates that a dclocalizcd two electron state 
is now possible in spite of the purely localized nature of 
the single particle spectrum and the repulsive character 
of the interaction. 

Of course, the most interesting question is to consider 
a system with a finite electronic density. As we have seen, 
because of the local Z2 symmetry, only pairs of electrons 
can propagate in a Z2 network. In particular, the single 
electron propagator — i(cj^(x, t)cj^(a;', t')) vanishes when 
i ji since it is not Z2 invariant, and it is very difficult 
to break spontaneously a local symmetrySi. This shows 
that the low-energy behavior of an electron fluid on a Z2 
network is not described by the Landau model of a Fermi 
liquid. As we shall see in section llVl below. it corresponds 
most likely to a new form of incoherent metal. But at this 
stage, the most natural thing to check is whether a super- 
conducting instability (induced by term [d] in Eq. 
occurs or not. 



B. Mean- field analysis of superconducting 
instability 

To begin this discussion, let us first consider a larger 
class of electron-electron interaction. The most general 
two-body potential interaction respecting the global spin 
rotation symmetry reads: 

^ = E EE /'"dx'y,^,?;.,(x,x') 

■ ■/ V / f -JO Jo 

4,o.,{x)4ro.'im,ro.'ii^')^^,o.,^^) (12) 

In this sum, {ij) and (i'j') are two pairs of nearest- 
neighbor nodes. As before, a (resp. a') labels one of 
the two hnks composing the loop {ij) (resp. {i'j'))- Let 
us use the cage basis, defined in Eq. l(Tn|l to rewrite this 
interaction. For this purpose, the first step is to express 
the local electronic density as: 

V'lar(2;)V'yar(^) = ^ (cL (a;)c„ (x) -f c]^ (a;)Cj.^ (x)) 

+ i(e-'^-°4(x)c^,(x) + e-'«--°4,(x)c,,(x)) 
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Because the flux through each elementary loop is half a 
flux quantum, e^'^'^° = — e~*^'3°. The interaction 112|) 
preserves the local Z2 symmetry if and only if the follow- 
ing conditions are satisfied: 

(x, x') = V^^, (x, x') = V^^, (x, x') (13) 

provided the pairs (if) and (i'j') are distinct, and: 

l^^^(x,x') = V^i;(x,x') (14) 

The first condition implies that 



when (if) and (i' f) are distinct, and the second condition 
that 



5^e-'^--l^^«;(x,x') = 

a, a' 



Therefore, all the terms which do not preserve the parity 
of the total electron number in each cage are eliminated. 
When both conditions (|13|1 . (|14|l on the interaction po- 
tential are satisfied, the interaction term H12|l now be- 
comes: 



EEE 



IJ,V 3 OLOL TT 



dx / dx' V^^,.,{x,x') 



cl{x)c^^(x) + cUx)c^^(x)] [4,,(x')c,,,,(x') + 4,,(x')c,,,,(x') 



(15) 



+ 5. 



e-'^-°4(x)c,,(x) + e-'^^'°ct (x)q,(x) e-'^'-'c,L(x')c,,(x) + e-''^^"' cUx)c„(x) 



We have used the standard notation : O : for the normal 
ordering of the operator O. 

Let us now consider the possibility of a superconduct- 
ing instability. In a variational approach, we would use 
a Hartree-Fock-Bogoliubov trial state. In the presence of 
pairing, the average kinetic energy increases, so it would 
have to be compensated by a lowering in the average 
interaction energy. As usual, this quantity involves aver- 
ages of the form {c^c ){c^c ) and of the form {c^c^){c c ). 
In the spirit of a mean field theory, we shall keep only 
the latter terms. Once again we also argue that the local 
Z2 symmetry is not broken, even in the hypothetical su- 
perconducting state. Physically, this is reasonable since 
the size of the Aharononv-Bohm cages being finite, there 
is a finite gap between the Z2 sector with = 1 and the 
sector with = — 1. As discussed below, this gap can be 
as large as a few Kelvin degrees for a GaAs Z2 network. 
With an unbroken Z2 symmetry, the superconducting or- 
der parameter (c|^(.t)c]^(x')) vanishes if i 7^ j. 

Let us first consider a singlet pairing. This leads to: 

{c\,(x)c\^(x')) = {clix)cl(x'))=0 
{cU^)cl(x')) = -{cl(x)cl^(x'))^X(x,x') 

In this case, the local contribution of Cooper pairing to 
the trial potential energy is: 

F,^,?;.,(x,x')|e-^'-A,(x,x') + e-5«^'°A,(x,x')P (16) 



In the case of triplet pairing, we have a complex vector 
Aj(x,x') defined by: 

AJ(x,x') = -L(-c],(x)c],(x')+4,(x)4,(x')) 
Ay(x,x') = -^{c],(x)cl^ix') + c],(x)c],(x')) 
A,^(x,x') = -^{c},ix)cl(x') + c],(x)c]^(x')) 

For triplet pairing, the orbital angular momentum of a 
Cooper pair is odd, so Aj(x, x') = — Aj(x', x). It is sim- 
ple to check that the above expression (|16|l still applies up 
to replacing the complex scalar Aj(x, x') by the complex 
vector Aj(x, x'). 

To conclude, we show that the part of the trial en- 
ergy which is directly sensitive to a Z2 invariant Cooper 
pairing amplitude is positive for a positive interaction 
potential compatible with the local Z2 symmetry. So, 
at least in a mean-field picture, a superconducting in- 
stability cannot occur. This statement was not a priori 
obvious, given the existence of delocalized states for the 
two-particle problem, where both particles are located in 
the same cage at any time, even with a repulsive inter- 
action. 
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C. Discussion for networks in semiconductor 
heterostructures 

For applications to real systems, it is important to 
present some order of magnitude estimates. For wires 
etched on a 2D electron gas at a GaAs/GaAlAs inter- 
face, with a typical carrier density of 3 10^^ cm~^ and a 
transverse width w = 100 nm, the number of transverse 
conduction channels N — k-pw/n is aroimd 4. For a regu- 
lar network of single channel wires, we have investigated 
interaction effects on the low-energy scattering matrix at 
the node92^. The renormalization of the transmission co- 
efficient r(A) as a function of a typical energy scale A is 
given by: 

where a ~ U{k ~^ 0)/{hv-p) is proportional to the Fourier 
transform of the interaction potential in the long wave- 
length limit. This result corresponds to nodes having 
a complete symmetry under any permutation of the in- 
cident wires. Note that it has the same form as for a 
single elastic scatterer in a Luttinger liquid^i^, or a 
single crossing connecting several semi-infinite Luttinger 
liquidsS^. For a single wire of width w and length L, 
U{k — > 0) ~ ^ hi{2L/w). Choosing typical values as 
TO* = 0.067 TOe and = 12.9 for GaAs, L = 2 jim and 
w = 100 nm, we obtain a = 0.4 which is a rather large 
value. This shows that the Coulomb energy is nearly as 
large as the kinetic energy of conduction electrons. 

As usual for ID electron systems with repulsive in- 
teractions, the leading instability occurs because of the 
divergence of the static charge and spin susceptibilities 
at the wave- vector q = 2kp and favors the formation of a 
spin-density wave. A derivation of this mean-field tran- 
sition on a Z2 network is sketched in Appendix 1X1 But 
with a = 0.4, the critical temperature for the onset of this 
mean-field instability is not much smaller than the Fermi 
energy (see for instance Eq. (|A1|) '). which shows that the 
system is still in a one-dimensional regime. Indeed, the 
cross-over from one to two dimensions takes place when 
the ID thermal length Lt = hv-p /k^T becomes compara- 
ble to the distance L between nearest-neighbor nodes. In 
other words, the corresponding energy scale k-^T should 
be as low as the energy level splittings inside Aharonov- 
Bohm cages, that is of the order of hv-p/L. Since Ap = 50 
nm {L/Xp = 40), the cross-over to 2D behavior appears 
at an energy scale one order of magnitude smaller than 
the mean-field Peierls instability scale. In practice how- 
ever, because of the small effective mass in GaAs, the low 
energy scale hvp/L can be of the order of a few Kelvin 
degrees. This scale sets also the order of magnitude of 
the Z2 gap already mentioned. 

It is well-known that in a purely ID system, fluctua- 
tions in the Peierls and in the Cooper channels have a 
tendency to neutralize each other, and the result is in 
general some kind of Luttinger liquid2&. As mean-field 



analysis fails to predict such a behavior, we should there- 
fore revisit the result of the previous section using an 
approach which is able to treat both Peierls and Cooper 
channels on an equal footing. This is exactly what the 
renormalization group (RG) is able to provide, at least in 
the weak interaction limit. To be completely honest, it is 
not at all clear that a pcrturbativc RG can give reliable 
quantitative estimates when a is as large as 0.4. But for 
the sake of establishing the presence or the absence of an 
instability, we expect it to produce a valid qualitative pic- 
ture, specially in the case, as we shall see below, where it 
predicts an instability to a gapped state already for small 
coupling. It might happen, as in some Hubbard models 
with extended internal symmetries'^, that some features 
of the gapped phase are modified when the coupling be- 
comes large compared to the electronic bandwidth. But 
in this example, the gap never closes when the system 
evolves from weak to strong coupling. 



IV. RENORMALIZATION GROUP ANALYSIS 
A. Low-energy model 

A complete RG analysis for our original system would 
be very complicated, because of the absence of transla- 
tion invariance inside Aharonov-Bohm cages. In partic- 
ular, network nodes induce Fricdel density oscillations 
in the electron fluid, which renormalize the associated 
single electron scattering matrices in the presence of in- 
teractions. This mechanism is precisely at the origin of 
the scaling shown in Eq. lfT7|l ^^i^^i^^. Physically, this in- 
teraction effect shows that wires have a tendency to be- 
come disconnected at low energies, which stops below the 
cross-over into the 2D regime. As shown previouslji^, an 
electronic wire network becomes an insulator for com- 
mensurate filling factors, or simply highly resistive, with 
a weakly dispersive band crossing the Fermi level for 
generic filling factors. However, for a system with local 
Z2 symmetry, this renormalization mostly modifies the 
internal dynamics inside a given Aharonov-Bohm cage, 
and not the general structure of the electron-electron in- 
teraction as it appears in Ea. 1111 At this point, we shall 
now assume that details of this internal dynamics inside 
Aharonov-Bohm cages do not play much role in the anal- 
ysis of interaction effects at low energy, provided the lo- 
cal Z2 symmetry manifested in Eq. (|ll|l is still respected. 
The validity of this assumption will be examined in more 
detail at the end of this section, once the results of the 
RG treatment are available. 

This leads us to consider the following model, where 
each Aharonov-Bohm cage is replaced by an infinite ID 
Luttinger liquid, described by the electronic fields c|^(a;), 
Cj^(a;), i being a cage label and t the spin projection. 
The internal energy levels of this collection of modified 
Aharonov-Bohm cages is naturally described by the fol- 
lowing kinetic Hamiltonian: 
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(18) 



where hvp has been set equal to unity. As usual, right 
and left moving fields have been introduced, labelled by R 
and L respectively, so that (a;) = ^ ^ (x) + j (x) . 
Note that replacing finite Aharonov-Bohm cages by in- 
finite ones is only valid in the ID regime, when we can 
neglect the discrete nature of the single particle spec- 
trum. The physical motivation for doing this is that as 
discussed in the previous section, the typical interaction 
scale is larger than the Z2 gap in present experiments. On 
the other hand, the perturbative RG used here requires 



small interactions in comparison to the single electron 
bandwidth. Because of the large value of L/Af, these 
two requirements are mutually compatible. 

Two nearby cages i and j are coupled by a two-body 
interaction subjected to the following constraints: it has 
to preserve the local Z2 symmetry or equivalently, to con- 
serve the parity of the total number of electrons in each 
cage, and to simplify the treatment, to preserve also the 
translational symmetry within cages. Both requirements 
are satisfied by the following Hubbard-like interaction: 



Hi, 



U / dx 



4] (2^)^^ {x)c^] {x) + c]^ {x)c]^ {x)c^^ {x)c^^ {x) + {x)c]^ {x)c.^ {x)c^^ {x) + c]^ {x)c\^ {x)c^^ {x)c.^ {x) 
+ 4T(^)sV^)"^a(^)S-T(^) + c]^(2;)cj^(a;)c^-^(x)c^^(x) + c\^{x)c\^{x)c^^{x)c^^{x) + c\^{x)c]^{x)c^^{x)c^^{x) 



Using the decomposition {x) = Cp , ^ (x) + cj^ j ^ [x) 
generates many terms, among which we shall neglect two 
subsets. The first subset contains all the purely forward 
scattering processes on a given side of the Fermi-surface, 
i.e. of the form cj^CpCpCpj or c^c^c^Cj^. In perturbation 
theory, they do not give singular terms in the low-energy 
limit, so we shall ignore them. The second subset con- 
tains all Umklapp terms, of the form CpCpCLCpj (resp. 
'^r'^r'^l'^l) which are relevant only for a filled (resp. half- 
filled) single electron band. The most interesting terms 
are of the type Cj^ClClCj^ which induce logarithmic di- 
vergences to the dressed two-particle interaction in the 



low energy limit. These terms are adequately treated by 
the renormalization group. Introducing now generalized 
charge and spin densities defined as: 

PR(L),„,«'(a;) = E4,(L),„,r(2^)^T,r'CR(L),„'.r'(2^)' (l^) 
r,r' 

S'R(L),n,n'(a;) =^4(L),„,r(2^)'^r,r'Cj^(L)_„,^^,(a;), (20) 
r,T' 

where cr = (cr"^, cr^, ct^) are the usual Pauh matrices, we 
may finally write the interaction term in the form: 



Hi, 



U 



/oo 
d.T^ : 2pR (x) — 2S'r 

-oo '- 



(2;)pL,n+i,«+i(a:) - 5'r,„.„S'l,„+i,«+i + PR,n+i,«+i(a;)pL,«.n(a;) - 

(2^)PL,n,n+l(a;) - 'S'r,„+1,„S'l,„,„+1 + PYl,n ,n+\{x) Pi. .n+\ ,n{x) - S'r 
(a;)pL,«+i,n(a:) - 'S'R,n+l,nS'L,n+l,n + PR,n,ri+l (2;)pL,n.n+l 

ix) - S 



(21) 



r 



Here, we have chosen a particular ID geometry, namely plied to any regular lattice of Luttinger liquid-like cages, 
a chain of loops. Although the RG approach could be ap- we shall from now on focus on this example. The re- 
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suits we shall present shortly suggest that the large scale 
lattice geometry does not affect much the physical prop- 
erties of such a system. 

B. Renormalization group equations 

We shall now briefly derive one-loop RG equations. We 
refer the reader to Ref. and references therein for gen- 
eral details about the RG technique. Here we focus on 
what happens in the space of cage positions, since this is 



what makes the originality of our work. The microscopic 
Hamiltonian I|18I21|I serves as the high energy initial con- 
dition of RG flows. The kinetic part of the Hamiltonian is 
unrenormalized at one-loop order, while the interaction 
part flows. During the renormalization process, the inter- 
actions will develop a non-trivial dependence in the cage 
positions. The local Z2-symmetry as well as the trans- 
lational symmetry of the initial Hamiltonian H18I21|I will 
remain symmetries of the renormalized Hamiltonian. As 
a consequence, the flowing Hamiltonian is constrained to 
be of the following form 



i^o(A) 

ifint(A) 



/ ^E ■ 1 E [^ji^)P^.nAl)Ph.'n+'r.n+ji-q) + ( A) ^R^n,™ (g)5L,„+7,„+T, (-?) 
+ C^{A)pji^n+-y,n{q)Pl^,n+f,n{-q) + (A)Sr^„+-,,„ (g)S'L,, 1+7,11 (-g) 



(22) 
(23) 



+D''{A)pn,n,n{q)Pi^,nA-q) + {A) S 

{q)SL,n,n{-q) 



Note that the hermiticity of the Hamiltonian implies 
that = C^^^\ and the Right/Left symmetry gives 



c(s) 



is) 



b: 



:(s) 



The action of the four 



types of locally Z2-symmetric couplings A^, Bj, and 
D (forgetting about their charge/spin nature) are repre- 
sented in Fig. [S] 

In Eqs. (|22I23|) we used the Fourier transform of the 
field operators which is defined by 



/OC 
da;cxp(-ifca;)cp(L),„,^(2;)- (24) 
-oo 



The Fourier transform of the generalized charge density, 
for example, is given by 



Pr(l) 



da;exp(i(7a;)pR(L)^„^„/(a;) (25) 

^ E 4(L),„.r(fc + '?)cR(L),„',r(fc)' 



All momenta have been restricted to be smaller in ab- 
solute value than the cut-off A, in order to regularize 



ultra-violet divergences. Note that the Fermi momen- 
tum has been set to zero, which is allowed since we do 
not consider Umklapp processes. We have not denoted 
any momentum dependence of the interactions, because 
in the RG treatment, all momenta are taken at the Fermi 
surface. 

The high-energy Hamiltonian is defined at scale Aq, 
and the initial conditions of the RG flow, at this scale, 
are found from (|21|l and read 



^^(Ao) 
^;(Ao) 
D\A^) 
D'{Ao) 



i?^(Ao) = q(Ao) = C/((5^,i + <57,_i) 
q(Ao) = -C/(<57,i + <57,_i) 



sliAo) 
2U 

-2U 



(26) 



The interactions are renormalized both by particle-hole 
and particle-particle contributions. Written with the di- 
mensionless scale / = ln(Ao/A), the RG equations are 



11 




FIG. 5: Schematic representation of the four interactions compatible with the local Z2-symmetry, given in Eq. Each 
dashed ellipse represents a cage. An empty (full) circle represents a destruction (creation) operator. Coupling A is an inter- 
cage interaction, coupling B a cage exchange, coupling C a pair hopping, and coupling D an intra-cage interaction. 
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(27) 
(28) 

(29) 

(30) 



In these flow equations, we have denoted 
Gi * G2 
Gi * G2 



Ph - ^[G'iG'^ + ^GlGl), 

- f GI G% + GI G% + 2GI G% 



ph 2 

for the particle-hole channel, and 



(31) 
(32) 

(33) 
(34) 



for the particle-particle channel. The RG flow equations 
can be simplified further, but we shall leave them in the 
form H27I30|I because it clearly shows the origin of all one- 
loop contributions. As an illustration, we show in Fig. 1^1 
the bubbles rcnormalizing the pair hopping term Cj. 
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Numerical solution of the Renormalization 
Group equations 

1. Strong coupling Hamiltonian 



The flow equations H27I30|I cannot be solved analyti- 
cally, but one can integrate them numerically with stan- 
dard numerical routines such as Runge-Kutta of order 
4. One nice feature of the local Z2-symmetry is that it 
constrains interactions in such a way that they can de- 
pend on at most one cage index. Thus the number of 
equations that have to be solved only grows linearly with 
the total number of cages N , allowing to consider large 
system sizes. 

With the initial conditions H2t)|l . wc found that the sys- 
tem goes to strong coupling. The flow diverges at a criti- 
cal RG scale Ic, indicating a phase transition at a critical 
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FIG. 6: Graphs renormalizing the pair hopping term C^. The left and middle graphs are particle-hole bubbles, and the right 
one is a particle-particle bubble. For the latter, (3 is summed over, and when /3 = or 7, the quantity Co that appears should 
be understood as D, so as to yield Eq. I|29^ . 



temperature of order of magnitude Tc = Aoexp(— /c). 
Apart from finite-size corrections, the numerical solution 
of the flow reaches the following fixed direction 



Al, = 0, 



B"^ ^ D'' ^ -iBl^ = -SD" > 0, 



CI 



(35) 
(36) 
(37) 



for any value of 7. The corresponding effective interac- 
tion Hamiltonian can be written (going back to real space 
and setting the strength of the interaction to a value g) 



eff 



/oo 
dx^ ; \ipK,n+i,n{x)p'L,n,n+i{x) 

(x). : (38) 



This Hamiltonian has an SU(7V) symmetry in the cage 
indices, and it corresponds to a particle-hole pairing in 
the triplet channel. Consequently the system undergoes 
a phase transition to a spin-density wave state, as was 
predicted by our mean-field analysis (see appendix^). 
This is more transparent on the equivalent form of (|38|) : 



HI 



off 



-2.9 



' m,n fji^' vv' 



(39) 



For an initial interaction \J = 0.1, this gives a ratio 
rc(RG)/Tc(RPA) ~ 10-3, which is very small. The RPA 
thus overestimates the critical temperature by three or- 
ders of magnitude for U = 0.1. For the typical value 
a = 0.4 {IJ = 27ra) estimated in GaAs networks, this 
reduction effect is not as dramatic, since here we ob- 
tain: rc(RG)/rc(RPA) ~ 0.76. This gives further sup- 
port to replacing finite Aharonov-Bohm cages by infinite 
ones in our simplified model, since the energy scale which 
emerges from the RG treatment is still in the ID regime, 
for a typical choice of parameters such as L/Ap = 40. 



2. Intermediate energy regime 

Another feature of the RPA is that the pair hopping 
terms C^^^^ remain nearest neighbor all along the flow, 
as can be inferred from Eq. H29() . when neglecting the 
particle-particle contribution. This is obviously not the 
case any more for the full RG equations. As the pair 
hopping terms are the only processes that can lead to 
a conducting behavior, it is natural to study what hap- 
pens to them during the flow, before they get completely 
suppressed by the spin-density wave. 

We found numerically that, except at the very begin- 
ning of the fiow, the pair hopping terms are very well 
fitted by the following law 



• (cR,m,^x(^)'''M,M"^L,m,Ai'(^)) ('^Ln,!^ (^)'''i',i^"^R,n,iy' (^)) ' C"^^^) (/) = (-l)'^r'^^^^ (Z) exp[ 



(41) 



since it is bilinear in the 2fcF Fourier component of the 
local spin density. A more detailed analysis of this phase, 
including a discussion of its collective modes can be found 
in Ref. iH 

The RG calculation however gives a much smaller value 
of the critical temperature than the mean-field calcula- 
tion. We have computed the critical temperature given 
by the full fiow equations (|27I30() , as well as for the RPA 
approximation of these, which consists in neglecting all 
particle-particle contributions. We found the following 
numerical result 



/c(RG) 



1.33 ^ , 0.64 
and ;c(R.PA) 



U 



U 



(40) 



This shows the incoherent nature of these pair hopping 
processes. As an illustration. Fig. |7|shows ln|C^(Z = 1)| 
as a function of 7, for N — 128 cages and interaction 
strength U = 0.1. We have performed such fits for many 
values of / so as to extract the scale dependence of the 
quantities r'=(")(/) and ^"^""^l) defined in Eq. 0. The 
results for the charge and spin couplings are nearly iden- 
tical, so that we only show the behavior of the charge 
quantities, see Figs. |S1 and 1^1 From Fig.|Hl it is clear that 
the pair hopping is suppressed during the fiow, but there 
is an intermediate energy regime (2 < Z < 12) where 
the decrease rate of the pair hopping is not too large. 
Furthermore, we see in Fig. |3 that the typical range of 
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FIG. 7: The numerical results of In |C^(/ — 1)| as a function 
of 7, for A'^ = 128 cages and initial interaction U = 0.1, are 
perfectly fitted by a straight line. 




FIG. 8: Evolution of the strength of the interaction Cj{l) as 
a function of the RG scale I, for A'' — 128 cages and initial 
interaction U = 0.1. The flow stops at the critical scale Ic 
denoted by a dashed line. 



the pair hopping increases during the flow, but it does 
not exceed one lattice spacing. The pair hopping thus 
remains local, contrarily to the cage exchange terms By 
which become of infinite range during the flow, as was 
explained in the previous subsection. 



D. Validity of the simplified model 

An important question arises: when can the Friedel 
density oscillations be ignored? In a real system, two 
physical processes are in competition: the wire discon- 
nection at network junctions due to the S'-matrix renor- 
malization, and the SDW instabiUty discussed earlier. In 
order to estimate the typical energy scale of wire discon- 
nection, we use Eq. (|17|l and take the point where both 
terms in the denominator are of the same order of mag- 
nitude. 



l/2a 



(42) 



The simplified model discussed in this section gives a 
faithful description of the initial Z2 network only if the 
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FIG. 9: Evolution of the range of the interaction Cj{l) as a 
function of the RG scale I, for A'^ = 128 cages and interac- 
tion strength U = 0.1. The flow stops at the critical scale Ic 
denoted by a dashed line. 



typical energy scale for the SDW instability is larger than 
Tc(S). From we have: 



fcBrc(RG) ~ Spcxp 



1.33 



This gives a rather precise criterion, namely: 

0.655 



Rq 1.33 

— < exp( ) 



(43) 



(44) 



Note that because of the one loop RG approximation 
used for both mechanisms (i.e scattering on Friedel oscil- 
lations, and SDW instability), this criterion is indepen- 
dent of the interaction strength, provided it is small com- 
pared to the electronic Fermi energy. We may say that it 
is of a purely geometric nature, since it mostly depends 
on the amount of electronic backscattering induced by 
the network nodes. For an illustration, let us consider 
the two types of node geometries depicted on Fig. 13 For 
single mode wires, we have found that the lowest values 
of the ratio Rq/Tq for perfectly symmetric X or Y junc- 
tions are equal respectively to 1 and 1 /4. So changing the 
junction type from X (Fig. El a)) to Y (Fig.Elb)) could 
affect deeply the nature of the low-energy state. 

In the case of a single ID quantum wire, the RG ap- 
proach gives a vanishing transition temperature and the 
spin density wave is replaced by a Luttinger liquid fixed 
point at low energy. For Z2 networks in the regime where 
Tc(RG) > T'c(S), inter-cage interaction processes destroy 
this Luttinger liquid, and because of the flow to strong 
coupling, most elementary excitations acquire a gap of 
the order Tc(RG). In the other regime, (Tc(RG) < 
Tc(S)), our simplifled model discarding Friedel oscilla- 
tions cannot capture the low energy physics below Tc(S). 
Below this scale, the network "breaks" into independent 
disconnected links, which are expected to behave as in- 
dependent finite segments of Luttinger liquids. 

How could we distinguish between these two situations 
in real experiments? It may not be very easy at a quali- 
tative level, just looking at transport measurements. Re- 



14 



sistivity is expected to rise sharply below the largest tem- 
perature scale between TclR-G) and T'c(S), but for very 
different reasons in each case. If Tc(RG) > T'c(S), the sin- 
gle electron spectrum is gapped, and charge transport is 
provided only via a single charge collective mode, which is 
likely to be quite sensitive to pinning by random impuri- 
ties. If rc(RG) < T'c(S), the disconnection of the network 
into independent segments clearly induces a very high re- 
sistivity. The difference between the two situations could 
be evidenced by non-linear dc-transport. Indeed, in the 
former situation, a pinned spin-density wave can be de- 
pinned by a large enough electric field, and should lead to 
a large non-linear dc conductivity above a finite thresh- 
old value for the driving field. In the later case, possible 
non-linear effects on the disconnection mechanism have 
not been investigated to our knowledge. But it seems un- 
likely that they would yield a sharp threshold as in the 
case of collective depinning. 

The above discussion shows that it would be crucial to 
estimate the ratio Rq/Tq for a given sample. In present 
experiments, although the number of channels Nch is low 
(about five), it is not exactly one. So we have a Z x Nch 
conduction matrix at each node connected to Z nearest 
neighbors. This clearly calls for an extension of the the- 
ory to the case of multichannel wires, which goes beyond 
the scope of the present work. But even for a strictly sin- 
gle channel network, a direct experimental determination 
of Rq/Tq may not be straightforward. In a dc transport 
measurement of an ideal system, the Landauer approach 
shows that the resistivity is dominated by the contacts 
between the network and the outer current leads. As 
shown in2S, Rq /Tq controls the dispersion of the mini en- 
ergy bands associated to the lattice superstructure. If 
Rq/Tq is small, these bands are weakly dispersive, and 
the effect of static random impurities, always present in 
a real system, is stronger than for a large value of Rq/Tq. 
But in general, we do not have an independent measure- 
ment of the disorder strength. So at this stage, the main 
prediction is that some Z2 networks should develop a 
SDW instability and some others break into indepen- 
dent disconnected segments at low energy. Although the 
mathematical criterion H44|) is quite precise and simple, 
it may be difficult in practise to decide if a given sample 
will fall in one or the other class without actually doing 
the experiment. 



V. CONCLUSION 

In this paper, wc have studied a class of networks (Z2 
networks) for which all single electron energy cigenstatcs 
are localized in Aharonov-Bohm cages when an exter- 
nal magnetic field corresponding to half a flux quantum 
per loop is applied. The particularity of these networks 
(in comparison to other geometries showing Aharonov- 
Bohm cages such as the dice lattice), is the presence of 
a large group of local Z2 symmetries, which correspond 
to reversing the electronic current in all the loops adja- 



cent to a given node. Interestingly, these symmetries are 
preserved in the presence of a large class of interaction 
potentials, including the Hubbard local interaction. The 
conserved quantum numbers associated to this symme- 
try are simply the parity of the total electron number in 
each Aharonov-Bohm cage. Because of this strong con- 
straint, only pairs of electron can tunnel from one cage 
to another. 

A first consequence, assuming that the local Z2 sym- 
metry cannot be spontaneously broken, is that such a 
system is never a Landau Fermi liquid, since the single 
electron propagator is short ranged in space. This estab- 
lishes a qualitative difference between Z2 networks and 
geometries such as the dice lattice. In the latter case, 
Aharonov-Bohm cages are no longer protected by a sym- 
metry in the interacting system, and therefore a Fermi 
liquid induced by interactions seems perfectly plausible. 
This extends the distinction established for bosonic sys- 
tems for which a single boson condensate is always de- 
stroyed on a Z2 network, whereas it survives on the dice 
lattice in spite of a large discrete ground-state degener- 
acy. 

Using a simple mean-field analysis, completed by an 
unbiased renormalization group (RG) approach, we first 
conclude that no superconducting instability is expected, 
in spite of the particle pair hopping term generated 
by electron-electron interactions preserving the local Z2 
symmetry. The RG analysis shows that long-ranged pair 
hoppings are generated at intermediate energies, but the 
typical hopping range remains of the order of the lattice 
spacing. This time, these results are in sharp contrast 
with the bosonic counterpart, where the absence of a sin- 
gle boson condensate leaves room for a perfectly stable 
boson pair condensateii, with very interesting proper- 
ties such as topological ground-state degeneraciesi^. The 
main result given by the RG is the presence of a strong 
coupling low-energy fixed point, dominated by the long- 
distance fluctuations of a local spin density wave order 
parameter. The appearance of this fixed point can be at- 
tributed to the presence of Aharonov-Bohm cages, since 
in the model considered here, independent links would 
behave as Luttinger liquids. 

We should note that the RG calculation has been done 
on a simplified model, where all the lattice points inside 
a given Aharonov-Bohm cage are equivalent. So we have 
not incorporated the physics of Friedel density oscilla- 
tion induced by the nodes, which are known to renormal- 
ize the single-electron scattering matrices at these nodes, 
as long as the system is in the one dimensional regime 
above the energy scale hv-p/L, which can be as large as 
a few Kelvin degrees in a GaAs/GaAlAs structure. The 
validity range of this simplified model has been formu- 
lated in terms of a very simple criterion on the value of 
the reflexion coefficient at the network nodes. It has to 
be low enough, otherwise, the disconnection induced by 
electronic scattering on Friedel density oscillations dom- 
inates over the tendency to build a SDW instability, and 
then, the low energy state appears to be a collection of 
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independent Luttinger liquid segments. 

This study raises some open questions regarding the 
transport properties of such incoherent Z2 metal. It 
would be very interesting to compute the temperature 
dependence of the conductivity. At temperatures above 
the spin density wave instability, it is likely to exhibit 
a power-law behavior, with a presently unknown expo- 
nent. In the strong coupling regime, a gap for single 
electronic excitations is expected, so most likely, the tem- 
perature dependence of the conductivity will show an 
activated behavior. Another fascinating question is the 
transport through hybrid structures such as a series of 
a normal metal, a Z2 network, and a normal metal. In- 
deed, for driving voltages below the Z2 gap (of the order 
of hvp/L), incoming (resp. outgoing) electrons will have 
a strong tendency to be injected (resp. emitted) in pairs, 
although the Z2 network is far from a superconducting 
state. Maybe the most direct way to demonstrate this 
strange phenomenon in such a poorly conducting system 
would be through shot noise measurements. 
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APPENDIX A: SPIN DENSITY WAVE 
INSTABILITY ON A Z2 NETWORK 

For the sake of simplicity, let us consider the Hubbard 
local interaction. For an infinite ID system, we decouple 



the interaction as: 

U ((^|(x)VT(x))^j(x)^^(x) + {i^lix)ij^ix))^p\ix)i;^ix)) 

Let us assume that {ipl{x)^^{x)) ~ TAcos{2kpx). In 
linear response theory, self-consistency determines the 
mean-field critical temperature T* from the require- 
ment J7x^°H2fcF,T*) = 1, where x'^^Hl^T) is the static 
charge susceptibility of the non-interacting electron 
gas at wave-vector q and temperature T. Because 
x(°)(2fcF,T) = {hvF)-'^\n{aEF/kBT), where a is a nu- 
merical factor of order unity, the mean-field temperature 
T* obeys a BCS relation: 

ksT* ^aEFcM-^) (Al) 

For a Z2 network, our starting point is Eq. (|11|1 . Again, 
we are looking for a Z2 invariant order parameter of the 
form: (V'iY(2;)V'iT(-^)) ~ '''Ai{x). With such an order pa- 
rameter, the quantum expectation value of Eq. Hll() taken 
for a Hartree-Fock state is equal to — U{Ai{x) + Aj (cc))^. 
We see that choosing a uniform spin density wave order 
parameter Ai{x) ~ A definitely lowers the system inter- 
action energy, and we expect a mean-field transition tem- 
perature of the same order of magnitude as for the single 
infinite wire with a repulsive local interaction. Note that 
this transition is uniquely driven by the first two terms 
in Eq. (|11|) . that is the intra and inter cage interaction. 
The last two terms, i.e. the cage exchange and the pair 
hopping processes do not play a role in stabilizing the 
spin density wave state. As seen in section Hvl in a more 
elaborate treatment, these last two terms have the effect 
to lower the critical temperature by comparison to the 
mean-field estimate. 
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